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We explore the dark matter and collider phenomenology of the minimal gauged U(l)g model, 
consisting of a leptophobic Zb gauge boson, and an accompanying Higgs Sb- By requirement of 
anomaly cancellation, the fermion sector naturally contains a dark matter candidate—a Majorana 
isosinglet x stabilized by an inherent Z 2 symmetry. The absence of evidence for Z prime dijet 
resonances at the LHC suggests that the scale of symmetry breaking is Ab > 500 GeV. Saturation 
of dark matter abundance together with limits on the direct detection cross section (dominated 
by Higgs exchange) constrains the Higgs mixing angle to |0| < 0.22. For small mixing angles of 
|0| < 10“®, the 0(10%) branching fractions of the fermion loop-mediated 5 'b —>■ 77, Z7, ZZ modes 
may provide clues about the fermion content of the model at the LHC. 


I. INTRODUCTION 


Classical equations of motion derived from the stan¬ 
dard model Lagrangian imply exact conservation of the 
baryonic and leptonic currents 

q i=e,fi,T 


In contrast with other conserved currents such as elec¬ 
tromagnetic, color, and weak neutral/charged currents, 
these currents do not couple to any known gauge fields. 
This observation has motivated a number of authors to 
speculate on the possibility that these currents sepa¬ 
rately source fields that may have so far escaped detection 
m- It is recognized, however, that these symmetries 
are quantum mechanically anomalous, reflected by the 
nonconservation of the corresponding currents 


= dnJi" = 

MB ML ^0^2 2 






which makes them unsuitable to serve as sources for new 
gauge fields. To alleviate this difficulty, the standard 
model baryonic and leptonic currents are each extended 
by introducing additional fermionic degrees of freedom 
with opposite anomalous contributions, such that the to¬ 
tals are anomaly free: 


tM _ 7^ j_ 

'^B,tot ~ ^B,ext 

_ 7^ j_ 7^ 

*^L,tot ^L,ext 

O /M „ O 7M _ q 


( 1 ) 


The form of the extra pieces Jg and Jg , that are 
built out of the new fermions, are model dependent. It 
should also be checked that the modified weak isospin 
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and hypercharge currents remain anomaly free upon the 
addition of new degrees of freedom 

^M'^SU(2) “ '^M'^U(l) ~ 0 ■ (2) 

The total currents Jeton '^Ltot can then be consistently 
coupled to gauge fields Zb and Zb, thereby enlarging the 
particle physics symmetry gauge group to 

SU(3)c 0 SU(2)^ 0 U(l)^ 0 U(I)b 0 U(1)l. (3) 

This line of thought has lead to interesting physical 
results such as an understanding of the stability of dark 
matter (without imposing an ad hoc discrete symmetry) 

S , and possible unification among all group factors in 
at scales much lower than the traditional grand uni¬ 
fication scale without rapid proton decay 

The earliest models satisfying Q and (^ contained 
charged fermions that significantly modified the nominal 
properties of the standard model Higgs boson, and were 
subsequently ruled out by recent experimental findings 
at the LHC (see [7] for a brief review). 

Currently, there are two closely related viable minimal 
models. They contain a new set of colorless fermionic 
degrees of freedom, called leptobaryons, carrying both 
baryon and lepton numbers as required by anomaly can¬ 
cellation. The model in |3] was constructed first, and 
its phenomenology was studied in [TllSl (see Appendix [A| 
for its formulation). In this model, the extended bary¬ 
onic and leptonic currents contain both vector and axial- 
vector parts 


= + (4) 

The axial charges are fixed by requirement of anomaly 
cancellation, but the overall normalization of vector 
charges is unrestricted and represent a free parameter 
of the theory. Throughout this paper, we refer to this 
extension as leptobaryon model VA after the form of its 
U( 1 )b,l couplings. 

Recently [5], a more systematic study revealed a sec¬ 
ond consistent model which contains fewer fermion multi- 
plets. In this model, the extended baryonic and leptonic 
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TABLE I. Field content of leptobaryon model A, with the 
fermions given in two-component left-handed notation. All 
fields are color singlets, and the hypercharge is normalized 
according to Q = -\-Y. Leptonic fields bracketed in paren¬ 
thesis Zl, S'l, and iZe./j.T are not considered for the present 
phenomenological study. 


currents have no vector part; so we call this extension 
leptobaryon model A. So far, there has been no study of 
its phenomenology. Our goal for the present paper is to 
initiate the exploration of dark matter and collider phe¬ 
nomenology. 

Throughout this article we only investigate the phe¬ 
nomenology of low scale U(l)g breaking, and neglect 
U(1 )l. This choice may be viewed within the leptobaryon 
model as taking the high-scale limit of lepton number 
breaking. However, there is no compelling reason not 
to consider 11(1)^ breaking near the electroweak scale; 
in fact, its phenomenology has been studied recently in 
the context of leptobaryon model VA liiin]. However, 
we leave the corresponding exploration in leptobaryon 
model A for the future. 

The rest of the paper is organized as follows. In Sec¬ 
tion [H] we formulate the model and analyze the ground 
state particle spectrum, and in Section [llT| we discuss phe¬ 
nomenological constraints on the model including vac¬ 
uum stability. In Section |IV[ we make a brief comment 
on baryogenesis in this model. The longer Sections 
and I VI I are on dark matter and collider phenomenology 
respectively. We draw conclusions in section [VHl 


II. FORMULATION OF LEPTOBARYON 
MODEL A 

A. Field content and Lagrangian 

The gauge-Higgs sector of the standard model is aug¬ 
mented by the addition of a U(l)g gauge boson Zg and 
a Higgs singlet 

C = - i sin(e)H^,Z^" 

+ + iD^SBrD>^SB - V{H, Sb) , (5) 


where sin(e) parametrizes the gauge kinetic mixing with 
the hypercharge gauge field B^. The covariant deriva¬ 
tives acting on the scalar fields are 

D^H = + igf- + \ig'B^)H , 

= {dfi -b , 

and the tree-level scalar potential is 

V{H, Sb) = - ^IS^^Sb 

+ biSI,SB)^ + aH^HSI,SB. (7) 

Four Weyl fermion multiplets, the leptobaryons, are in¬ 
troduced and their charge assignments (see Table [I]) are 
fixed by anomaly cancellation. The scalar field Sb is 
assigned a baryonic charge of -1-3 to allow for Yukawa- 
type couplings with the leptobaryons, which in two- 
component spinor notation are given as 

C = -y^SI,^^-^SBf^-f^-^SBXX + c.c., ( 8 ) 

and would generate their masses through the Higgs mech¬ 
anism. In general, the standard model Higgs doublet can 
also couple to the leptobaryons 

C = -Xi{H^)x-X2iH^^)x 

- A3i7l'(f-S)T'-A 4 ^(f-E)iJ-bc.c.. (9) 

These interactions mediate the decay of the heavier lep¬ 
tobaryons into lighter ones via the emission of Higgs 
quanta. 

It is useful to note that the pure Sb Yukawa interac¬ 
tions in ([^ exhibit reflection symmetries corresponding 
to the transformations 



vk and ¥ -> -¥ 


to 

E ^ -S 

(10) 

(^2)x : 

X ^ -X- 



The addition of standard model Higgs Yukawa interac¬ 
tions in ([^ explicitly breaks this symmetry leaving a sin¬ 
gle accidental Z 2 symmetry corresponding to the com¬ 
bined reflection of all leptobaryons, and preserves the 
stability of the lightest leptobaryon. These interactions 
contribute to custodial isospin violation, and can also 
provide a source of CP violation. 

B. Electroweak/baryonic vacuum structure and 
particle spectrum 

Phenomenological viability of the model requires that 
both the standard model Higgs doublet and the baryonic 
Higgs singlet acquire vacuum expectation values (vEVs) 
in the ground state 

+ ^' + 

Sb = :^{vb + s' + #b) , 


( 11 ) 
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FIG. 1. An illustration of the mass spectrum of leptobaryon 
model A considered in this study. We take the isosinglet x Si'S 
the lightest leptobaryon, which is automatically stable by the 
Z 2 symmetry in (lOl. 


where v = 246 GeV, and wb is the VEV of the baryonic 
Higgs conventionally taken to be positive by the Z 2 sym¬ 
metry of the potential. 

The Higgs portal coupling —aWHS^SB induces mix¬ 
ing between the neutral components of the doublet and 
singlet Higgs fields. Our convention (opposite to that in 
0 ) for the transformation to mass eigenstates is 



/ cos 9 sind\ /h\ 
\^— sin0 COS0/ ’ 


-f <0< 


TT 

4 • 


( 12 ) 


It is a straightforward exercise to apply the minimiza¬ 
tion and local stability conditions to exchange the five 
potential parameters A, b, a in terms of the VEVs 

V and vb, masses mn and ms, and the mixing angle 
9 (see Appendix E- Global tree-level vacuum stability 
(boundedness from below) requires that A > 0, & > 0 
and 4A6 — > 0, which are automatically satisfied by 

taking positive VEVs and masses. 

Upon spontaneous symmetry breaking of 17 ( 1 ) b , the 
baryonic gauge field Zb acquires a mass given by 

'mzB = ^9 bVb , ( 13 ) 

and the Yukawa interactions in Q lead to masses of the 
leptobaryons 41, E and x given by 




V2 


VB . 


(14) 


The Weyl degrees of freedom are organized such that 
when expressed in terms of eigenstates of the charge op¬ 
erator, both members of the isodoublet 41 = and 

the charged component of the isotriplet S+ are Dirac 


fermions, and the neutral member of the isotriplet E*^ 
and the isosinglet x Majorana fermions. Note that 
the leptobaryons carry the same electroweak quantum 
numbers as the gauginos and Higgsinos in the minimal 
supersymmetric standard model. 

In what follows, we take the Majorana fermion x to 
be the lightest leptobaryon (stabilized by an inherent 
Z 2 symmetry) and therefore as a dark matter candidate. 
Electroweak self-energies induce mass splittings between 
the charred and neutral components of the isomultiplets 
given bj0 

Tn^-\- TTl^O 

= « 341 MeV , , 

2 (15) 

TOy;± — m^o 

= I [a mz — aw(wz — mw)] ~ 166 MeV. 

This allows the charged leptobaryons to decay to their 
respective neutral partners via a charged pion emission. 
The resulting particle spectrum is sketched in Fig. [I] 
Finally, the Higgs interactions with the leptobaryons 
in 0 lead to mixing among the multiplets, and mediate 
the decay of all leptobaryons to the lightest one via Higgs 
boson emission. We shall simplify the phenomenology by 
taking the couplings Ai... A 4 to be smalQ even though 
this is not phenomenologically necessary. 


C. Electroweak and extended baryonic cnrrents 

After investigating the scalar interactions of the lepto¬ 
baryons we now turn to the gauge interactions. In our 
notation, the gauge interaction Lagrangian of the lepto¬ 
baryons is given by 

— Jqq + C.C.) — gBZBfiJB I 

where in four-component notation, the extended elec¬ 
troweak currents are 

JNc = il- 

Jcc = ^^7'^^° - ^7'^SO , ( 16 ) 

and the extended baryonic (axial-vector) current is 

•^B.ext = -i (^7'^75^+ + rrisr) 

+ I (^7^75E+ + i^7'^75E0) + |x7'^75A • (17) 


^ Here, we cite the values given in m- 

^ Note that this assumption is compatible with technical natural¬ 
ness because the symmetry in is recovered in the limit where 
these couplings vanish. 
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The couplings of leptobaryons to electroweak gauge fields 
radiatively influence the production and decay of the 
new scalar state ■ In section |VI B[ we illustrate how 
this can make it possible to distinguish model A from 
model VA. 


III. CONSTRAINTS ON LEPTOBARYON 
PARAMETERS 


We shall take the kinetic mixing to zero for simplicity 
(consistent at the 2a level). 

Finally, the observation of a Higgs-like resonance by 
the ATLAS [T3] and CMS [15] experiments consistent 
with signal strengths predicted by the standard model 
puts a limit on any mixing that the standard model Higgs 
might have with other degrees of freedom. We adopt 
the 95% confidence limit on the mixing obtained by the 
ATLAS Collaboration [T5] 


Leptobaryon model A brings in a number of new pa¬ 
rameters beyond the standard model: the gauge cou¬ 
pling, mass and kinetic mixing parameter 

ttB = mzB, sine; (18) 

47r 

the baryonic Higgs mass and its mixing with the standard 
model Higgs 

ms, 0; (19) 

and the leptobaryon masses along with the Yukawa cou¬ 
plings to the standard model Higgs 

TO^, ruE, m^, Ai,A2,A3, A4. (20) 

In this section, we review some important experimen¬ 
tal and theoretical constraints that will motivate certain 
simplifying assumptions on these parameters leading to 
simpler phenomenology in the sections to follow. 

A. Electroweak precision 


|6i|< 0.35. (23) 

B. Vacuum metastability at one loop 

Similar to the destabilization of the standard model 
Higgs potential by the large top quark mass the 
leptobaryons are also expected to destabilize the poten¬ 
tial along the baryonic Higgs field direction, yielding a 
constraint on the mass spectrum. A complete vacuum 
metastability analysis involving the full one-loop correc¬ 
tions to the scalar potential is beyond the scope of this 
paper. However, we can find fairly robust constraints 
by studying the behavior of the potential in the leading 
logarithmic approximation. 

We solve the one-loop MS renormalization group equa¬ 
tion for the baryonic Higgs quartic self-coupling 

= ^[- 44 -32/1;-2/^+ 20^+4865^ 

+ 2&(-544 + 44 + 32/1 + vD + 20b^] , (24) 


Because the couplings of leptobaryons to the elec¬ 
troweak gauge bosons in Eq. (16) are pure vector, they 


nominally make no oblique corrections at one loop {S — 
SsM = T — Tsm = 0). However mixing among the lep¬ 
tobaryons induced by the Yukawa couplings Ai... A4 in 
Eq. @ leads to deviations suppressed by the baryonic 
Higgs and leptobaryon masses [T^. As mentioned near 
the end of Section |HB[ we take these Yukawa couplings 
infinitesimally small only to permit the decay of the heav¬ 
ier leptobaryons to the lightest one. 

Next, we consider the effect of kinetic mixing param¬ 
eter sin(e) on low energy physics. For the range of 
masses we consider in this study (~0.5—2.5 TeV), kinetic 
mixing is most constrained by its effect on the relative 
strengths of the weak charged and neutral currents [15] 
parametrized by the p-parameter. In the limit of small 
kinetic coupling. 


Ap = 


4. 


tan^(6>w)e^-f 0(^4^) + C>(e^) . (21) 


m 




Fits of the p-parameter to electroweak precision observ¬ 
ables suggest a deviation from unity po = 1.00040 ± 
0.00024 yielding the preferred range of 


0.29 < 




-TeV < 0.58. 


( 22 ) 


with the initial conditions given at po = ijns + 'rnz^)l2., 
the defining scale for the masses and couplings. The scale 
p = Adestab. at which the quartic coupling crosses zero 
indicates the onset of destabilization in the S'b direction, 
and represents the scale at which new physics must enter 
to restore stability of the scalar potential. 

We impose self-consistency on the model parameters 
by requiring that the potential be stable at least up to 
the largest mass in the theory, the leptobaryon masses 
We assume degeneracy among the leptobaryons 
so that //p = 2/E = neglect the running of all 

other couplings in p4] ) which is justified a posteriori due 
to the very low scale at which b crosses zero. 

By retaining the dominant Yukawa coupling contribu¬ 
tions (first three terms) in the right-hand side of ( [24| , 
the renormalization group equation can be readily solved 
in closed form yielding an approximate constraint on the 
particle masses 

< 0.86(^=«)''f (25) 

^ 9b ^ 

This approximation agrees with the numerical determi¬ 
nation to better than 5% and indicates that unless the 
gauge coupling pB is exceptionally small, fermion masses 
should be of the same order of magnitude as the baryonic 
Higgs and gauge boson masses. In particular, vacuum 
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stability imposes a stringent upper limit on dark matter 
mass. 


C. On the naturalness of small parameters 

In the current study, we take the kinetic mixing pa¬ 
rameter sin(e) to vanish, and (for the collider study) we 
consider the scalar mixing angle 9 reaching ^ 10“® in 
magnitude. The smallness of these parameters is not as¬ 
sociated with the emergence of any symmetry, and in this 
section we ask whether these choices are consistent with 
the radiative corrections. 

The kinetic mixing angle sin(e) is the coefficient of 
a marginal operator , and is therefore logarith¬ 

mically sensitive to the UV scale. Without informa¬ 
tion about UV physics, no meaningful statements can be 
made about the size of radiative corrections. However, if 
(setting aside the issue of vacuum stability) grand uni¬ 
fication as suggested in [5] is assumed and the U(l) hy¬ 
percharge and baryon factors are orthogonal, sin(e) can 
be estimated based on the renormalization group evolu¬ 
tion equation starting from sin(e) = 0 at the unification 
scalej^ Agut = 1-24 x 10^^ GeV. The quarks and charged 
leptobaryons (^/> and E) contribute to the beta function 
at one-loop order, given by 

+ (f + f 5 |) sin e] . (26) 

By evolving the kinetic mixing parameter down to the 
leptobaryon masses (^TeV scale), we find sin(e) « 0.06. 
We emphasize that this estimate should be taken with 
a grain of salt given that vacuum instability sets in at a 
much lower scale. New physics that stabilizes the poten¬ 
tial is likely to modify the beta function. 

The scalar sector mixing angle 9 is derived from the 
Higgs portal coupling a, which is also the coefficient of a 
marginal operator HS^Sb- In principle, a renormal¬ 
ization group evolution from the unification scale would 
provide an estimate at the TeV scale. But in this in¬ 
stance, more information about the structure of unifica¬ 
tion in the scalar sector is needed to make such an esti¬ 
mate. However, we mention that within the low energy 
leptobaryon model, renormalization effects are propor¬ 
tional to itself at the one-loop level. Independent renor¬ 
malization starts at the two-loop level, given by 



Therefore, we expect that the mixing angle is fairly stable 
against quantum corrections. 


IV. BARYOGENESIS 


In O [7], it was shown that if U(I)l is broken at a 
much higher scale, the decay of right-handed Majorana 
neutrinos can create an initial B — L asymmetry in the 
early Universe. The leptobaryons tp and S carry weak 
isospin, and therefore modify the ’t Hooft operator [5] 
associated with the sphaleron 

. (27) 


Above electroweak and baryon symmetry breaking scales, 
the modified sphalerons process the initial asymmetry 
A{B — L)sm into a final baryon asymmetry ABj^, with 
a conversion factor given by 

32 

= dd^(^ - -^)SM « 0.32 A{B - L)sm ■ (28) 

■’ yy 

This is slightly different from the standard model conver¬ 
sion factor first computed in m 

OQ 

= ^A(H - L)sm « 0.35 A{B - L)sm ■ (29) 


Provided that H acquires a VEV earlier than (or at the 
same time as) 5 'b does in the history of the Universe, 
the sphaleron processes decouple, and the relationship 
in (281 is frozen in. In view of the hierarchy of scales 
assumed in this model, it seems more plausible that H 
and acquire VEVs simultaneously. 

We mention here an important caveaf]^ that was not 
fully appreciated in [S]. If ^b acquires a VEV before H 
does, the Dirac mass term for imposes the chemical 
equilibrium condition = /i^. Then, the conversion 
factor for the final baryon asymmetry reverts back to the 
standard model one in (29). 


V. LEPTOBARYONIC COLD DARK MATTER 


We begin our phenomenological exploration of lepto¬ 
baryon model A by investigating the properties of the 
dark matter candidate. There are three electrically neu¬ 
tral leptobaryons, and y, each one of which could 

in principle account for the observed dark matter abun¬ 
dance. In this paper, we take the isosinglet y as the 
lightest stable leptobaryon to serve as the dark matter 
candidate, and defer the exploration of the isotopically 
charged leptobaryons ip° and E° as dark matter to future 
study. 

Note that since y is a Majorana fermion, its coupling 
to Zb is purely axial [see 0] , and is in contrast with the 
dark matter candidate of model VA which also couples 


vectorially [see (A2)]. This leads to major qualitative dif¬ 
ferences in the dark matter phenomenology with respect 
to model VA, largely because several possible scattering 
channels are velocity suppressed. 


® See Table I of [6] 


^ We thank David Morrissey for alerting us to this possibility. 










A. Dark matter relic abundance and direct 
detection 
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We assume that the dark matter is produced by stan¬ 
dard thermal freeze-out. The dark matter abundance in 
the Universe is governed by the Lee-Weinberg equation 

m 


AY 

— =Z{x)[Y,l{x)-Y^x)] 


(30) 


where Y(x) is the dark matter number density normal¬ 
ized by the entropy density, and x = m^/T is the scaled 
plasma temperature. The factor Z{x) is given by 


/ TT m^TOp 

= y ’ 


where the thermally averaged cross section is 

(va) = f 
Ji 


y/s{s - Am\)a{s)Ki{Y^)x 
4 m, ?>m\Kl{x) 


(31) 


(32) 


The large x limit of the solution to ( |30[ ) determines the 
thermal relic abundance. 

We numerically integrate the thermally averaged cross 
section, retaining all orders in the velocity expansion to 
accurately capture the resonant and threshold behaviors, 
and we solve (30) in the freeze-out approximation |21| to 
give 


^(o®) = [yY) + ■1'^ ^ ^ 


with the freeze-out temperature x / determined by match¬ 
ing. To improve computational speed, we neglect the first 
term in (331 which leads to an overestimation of the relic 


abundance by no more than 5%. The thermal relic abun¬ 
dance is then obtained from 


Udm = "^r(oo), 


(34) 


where sq = 2970/cm^ is the entropy density today and 
Pc = 1.05394 X 10~^h^ GeV/cm^ is the critical density. 

In Table [nj we list the tree-level dark matter annihila¬ 
tion channels, all of which are taken to compute the dark 
matter relic density. Channels going directly to standard 
model particles such as xx HSb and XX ^ HH are 
omitted because they are suppressed by at least 1 order of 
magnitude owing to the collider constraint on the scalar 
mixing angle |0| < 0.35. To facilitate the identification 
of velocity suppressed amplitudes, the lowest contribut¬ 
ing partial waves are also given. Amplitudes for channels 
with an initial P (or higher) wave will be down by at 
least one power of dark matter velocity. 

In Fig. we plot the dark matter relic abundance 
as a function of dark matter mass for a 
typical choice of model parameters. The observed 
abundance at = 0.12 [52] is indicated by the 


Resonant annihilation channels 
XX m 


^P+ ^1- ^ qq 

Non-resonant annihilation channels 


XX ->■ ^bSb 



TABLE II. Set of dark matter annihilation channels in the 
limit of vanishing scalar H-Sb mixing angle 9, with the low¬ 
est nonvanishing partial waves indicated. Channels with an 
entering S-wave (xx ZbSb and xX ^b.^b) are not ve¬ 
locity suppressed. 



FIG. 2. Dark matter relic abundance (black) as a function 
of the dark matter mass for hxed Ob , , ms and 

9. The red line corresponds to the measured abundance 
= 0.12. The regions labeled ‘resonant’ and ‘nonreso¬ 
nant’ indicate the principal dark matter annihilation channels 
(see Table 0 that determine the relic abundance. 


red line, and intersections determine values of that 
saturate the relic density. The two intersections to the 
very left correspond to comparatively light dark matter 
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FIG. 3. Processes contributing to the spin-independent direct 
detection cross section, with the lowest nonvanishing partial 
waves indicated. Scattering mediated by exchange is ve¬ 
locity suppressed because it is a (parity-change) A£ = 1 tran¬ 
sition, while scattering mediated by H and 5 b exchanges is 
mixing angle suppressed. 


— 11-12; 0 —>0 



(m^ « ttizbI'^), and are determined by the resonant 
annihilation channel XX QQ- The single inter¬ 

section to the right gives the correct relic abundance for 
heavier dark matter masses, and arises from nonresonant 
annihilation which is usually dominated by the velocity 
unsuppressed channels xX ■^b5'b and xx —>■ ZbZb- 
To simplify the phenomenology, we will later separate 
the discussion corresponding to the two distinct mass 
ranges of the dark matter candidate. 


At tree-level, the dark matter-nucleon scattering chan¬ 
nels for direct detection are the t-channel Zb, H and 5 b 
exchanges. The Feynman graphs are shown in Fig. [^to¬ 
gether with their lowest contributing partial waves. The 
Zb exchange is velocity suppressed since it involves a 
change in the orbital angular momentum on account of 
the axial coupling to dark matter. On the other hand, the 
H and 5 b exchange is suppressed by the mixing angle. 

At leading order in the velocity expansion, the lab 
frame spin-independent direct detection cross sections, 
in the limit ^ mj\j, are 


o-SiI^b) =187ra|'^^u^ , 

/c ^ 72GFaB .2a 2a 

o'si(5b) =- -j= —sin t^cos u 

V2 

mlm%( 1 

^ 




(35) 


(36) 


where ruAr is the nucleon mass and = 0.32 is the 
effective nucleon matrix element taken from lattice cal¬ 
culations [53]. The total direct detection cross section is 
given by the sum ctsi = o'si(^b) + crsi(5B). Orthogonal¬ 
ity of the partial waves forbids an interference term. 

Altogether, the dark matter phenomenology princi¬ 
pally depends on five free parameters 


mzs, rns, m^, ub, 0. (37) 


B. Resonant dark matter 


Motivated by the vacuum stability bound (25), we 


start by assessing the prospects for the relic abundance 


FIG. 4. Gontours of fixed mz^ (brown lines, in GeV) and 
Ob (red lines) in the plane to achieve dark matter 

relic abundance saturation, for in the resonant annihila¬ 
tion region of Fig. ^ (blue dot). In the upper-right region, a 
nonperturbatively large coupling is needed. The limit 0 > 0 

is taken to decouple the scalar exchanges from direct detec¬ 
tion. The current LUX (2013) bound [2] and the XENONIT 
projections m are shown. 


being saturated by light dark matter due to the resonant 
(velocity suppressed) annihilation channej^xX Zg ^ 
qq in the early Universe. The annihilation cross section 
is given by 


tr = 97raB\/1 — 


mzj,Tzj,BT{ZB X) gg) 

s (s - 


(38) 


where th e Zb total width F Za is evaluated at s = 'rn^_^ 
given in (Cl). The annihilation cross section depends 


only on ub, wzb m^. For the moment, we decouple 
the scalar particles H and 5 b from the problem by taking 
0 —)■ 0 so that the only direct detection channel is via the 
t-channel Zb exchang^ We return to discuss the effects 
of restoring the scalar exchanges below. 

In this limit dependence on ms and 9 in the x~^ cross 
section drops out, and a positive signal by direct detec¬ 
tion experiments can fully determine the three remaining 
parameters in (37). For a given dark matter mass 


the direct detection cross section in (35) is inverted to fix 


ttB as a function of . The gauge boson mass mz^ is 
then determined by requiring saturation of the observed 
relic density. For illustrative purposes, we take the lower 
dark matter mass as shown in Fig. [^by a blue dot. 

In Fig. [^ we display the contours of constant cou¬ 
pling strength q;b (in red) and gauge boson mass mz^ 


® The resonant annihilation channel xX — ^ H/Sb — > QQ is both 
velocity and mixing angle suppressed. 

® Direct detection experiments report a limit on trsi under the as¬ 
sumption that to leading order it is velocity-independent. Since 
the Zb exchange gives a cross section that is velocity-dependent 
to leading order, we scale the (tsi{Zb) in (|35| l by the isothermal 
halo model of dark matter distribution 1261 to obtain an effec¬ 
tive cross section cto ~ (7p{u}/(u“^}. The net effect is to use an 
effective lab-frame velocity of u 0.0093c. 
















(in brown) in the m^-asi plane assuming relic abundance 
saturation. The shape of the brown mz^ contours reflects 
the shape of the relic abundance curve in Fig. [^in the res¬ 
onant region. The upper-right area (large and crsi) is 
where a nonperturbatively large ae is needed to achieve 
relic abundance saturation. The strongest direct detec¬ 
tion cross section bound by LUX (2013) [Ml currently 
excludes comparatively light dark matter my ^ 150 GeV 
assuming perturbative couplings and 9 = 0. Projections 
by the XENONIT Collaboration [25] indicate sensitiv¬ 
ity to larger dark matter masses reaching ^ 500 
GeV. We find that for nonvanishing 9, the velocity unsup¬ 
pressed H and Sb exchanges (see Fig. contribute to 
direct detection, but do not lead to significantly stronger 
exclusions. 



FIG. 5. The spin-independent direct detection cross section 
asi (dashed blue) as a function of the dark matter mass 
for several values of the Higgs mixing angle 9, for fixed mz^, 
ms and ob- The LUX (2013) bound is shown in black. At 
the dark matter mass mx determined by relic density satu¬ 
ration (red line), the upper limit on 6 is obtained when ctdd 
intersects the LUX bound. 


C. Nonresonant dark matter 

We now turn to the case of heavy dark matter (cor¬ 
responding to the nonresonant region in Fig. where 
rrix ^ We find that a much larger region in 

parameter space is allowed, and that current experimen¬ 
tal bounds on the spin-independent direct detection cross 
section provide limits on the mixing angle 9 that are more 
stringent than from the Higgs signal strength measure¬ 
ments at the LHC, which we discuss here. 

The annihilation channels involving the baryonic Higgs 
are open in this scenario, and play a crucial role in deter¬ 
mining the relic abundance. Thus, the dark matter phe¬ 
nomenology is influenced by all five parameters in (371. 


At leading order in the velocity expansion, the thermally 
averaged nonresonant annihilation cross sections are 

{vaixx^ ZbSb)) = ^ ml 


64 mtm% 

X 


{vcrixX ZbZb)) = 


SIttqI (m^ - 

{rn^z^ — 2m^)^ ’ 


(39) 

(40) 


where A(a, b, c) = + b'^ + — 2ab — 2ac — 26c. 

In view of the LHC upper bound on 0, we begin by 
taking 0 = 0 in the annihilation formulas (this approxi¬ 
mation will later be shown to be very good). Then, for 
fixed values of ms and obj we determine the dark 
matter mass by requiring saturation of the relic abun¬ 
dance. These values are shown as solid black contours in 
Fig. [^with ms = 1 TeV (left) and a. = 0.005 (right). 

The direct detection scattering cross section is now 
dominated by the t-channel H and exchanges given 
by Eq. (36), with the velocity suppressed Zb contribu¬ 


tion being negligibly small. The dependence of the cross 
section on the scalar mixing angle is particularly impor¬ 
tant. In Fig. the behavior of the spin-independent 
cross section (in dashed blue) is sketched for several val¬ 
ues of 9 along with the LUX (2013) [Mj upper limit in 
black. For the value of m^ determined by relic abun¬ 
dance saturation (red vertical line), direct detection null 


results are used to find an upper bound on the allowed 
mixing angle, determined by the intersection (indicated 
by a red dot). 

The contours of constant upper limits on the Higgs 
mixing angle from the LUX (2013) data are shown as 
dashed blue lines in Fig. [^with labels outside parenthe¬ 
sis. The weakest upper limit is around 9 = 0.22, and 
justifies a posteriori the approximation 0^0 originally 
taken in the annihilation formulas. This procedure can be 
viewed as the first step of a rapidly convergent iterative 
process to accurately determine the upper bound on 0 set 
by direct detection experiments. Shown in parenthesis on 
the same contours are the expected limits on the mixing 
angle derived from the XENONIT projections, and are 
significantly stronger than the LHG bound (|0| < 0.35). 

The rise in dark matter mass m^ required for satura¬ 
tion as the coupling an is lowered (left panel of Eig. 
can be qualitatively understood from Fig. A decrease 
in Ob implies smaller annihilation cross section, which 
leads to larger relic density corresponding to an overall 
upward shift in the black curve in Fig. This then 
displaces the intersection point (green dot) to the right 
implying heavier dark matter mass. The rise in m^ with 
increasing mz^ and ms (right panel of Fig. is more 
readily understood as requiring the nonresonant annihi¬ 
lation channels to remain open. 

With such heavy dark matter x even heavier 
charged leptobaryons ^ and E, the bounds derived from 
one-loop vacuum stability in Section HI B| become impor¬ 
tant. Displayed in the same plots is the boundary (red 
line) where the vacuum becomes unstable at a scale equal 
to the dark matter mass. In the region to the right of the 
red line, the vacuum becomes unstable at a scale below 
the dark matter mass, and signals the presence of new 
physics to restore vacuum stability. 

Despite bounds from one-loop vacuum metastability, 
the allowed region of parameter space for larger dark 
matter masses m^ ^ ’’mzs — is much wider than 
that for lighter dark matter masses ^ /2 ■ In ad¬ 

dition, direct detection experiments are becoming more 
sensitive to the Higgs mixing angle and we expect that 
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ms = 1000 GeV «B = 0.005 




FIG. 6. Contours of constant dark matter mass (black, in GeV) saturating the thermal relic abundance in the ‘nonresonant’ 
annihilation region of Fig. (green dot) with the approximation 6 = 0. Dashed blue contours indicate upper limits on the 
Higgs mixing angle 9 obtained from direct detection spin-independent cross section bounds. Values outside parenthesis indicate 
current limits based on the LUX (2013) data [24], and the values within parenthesis are the forecasts based on XENONIT 
projections |^. The red line demarcates the region below which the vacuum is radiatively unstable at a scale lower than the 
dark matter mass (see Section IIIB I. Left panel ms = 1000 GeV fixed, and right panel ctb = 0.005 fixed. 


future experiments will be able to tighten the limit on 
6 significantly more strongly than that from the Higgs 
signal strength measurements at the LHC. 


D. Upper bound on U(l)g breaking scale? 

We close this section by revisiting the arguments in [5| , 
in the context of leptobaryon model VA, leading to an 
absolute upper bound on the baryon symmetry breaking 
scale (^ 35 TeV) that is much lower than the one derived 
from unitarity (^ 300 TeV) |17|. Such a bound would 
be welcome, as it would imply falsifiability in the near 
future. In an attempt to derive a similar bound in the 
present model, we found that the arguments are based on 
an invalid implicit assumption, and upon further scrutiny 
the upper limit is lost. 

The implicit assumption made in the argument is that 
among the possible dark matter annihilation channels, 
the resonant channel XX qq is always the most 

efficient. It is then shown that as the symmetry breaking 
scale, i.e. the Zb mass, is parametrically increased, the 
thermal annihilation cross section for this channel drops. 
Eventually, the magnitude of the cross section becomes 
too low, unavoidably leading to an overdosed Universe. 
The value of when this occurs evidently becomes 
the upper limit on the baryon symmetry breaking scale. 

However, we point out that even for large mzg, the 
nonresonant annihilation cross sections can be made ar¬ 
bitrarily large by increasing ob thereby guaranteeing suf¬ 


ficient dark matter annihilation. One may worry that a 
nonperturbatively large ob or is needed, or that uni¬ 
tarity must be violated before reaching sufficient annihi¬ 
lation. However, this concern has already been addressed 

in ETj. 

The arguments above also rely on y being the lightest 
leptobaryon. It is possible for the isodoublet ijj to be the 
lightest leptobaryon, with the electrically neutral com¬ 
ponent as the dark matter candidate. Because of its 
electroweak charges, many additional annihilation chan¬ 
nels are open, and it is not difficult to imagine a scenario 
with high scale U(l)g breaking and where '0® forms a 
subdominant part of multicomponent dark matter. 

We mention that, within the minimal leptobaryon 
models, it may be possible to use bounds from vacuum 
metastability in conjunction with an analysis of the anni¬ 
hilation channels to place an upper bound on the U(l)g 
breaking scale. However, given that new physics is ex¬ 
pected to enter, the bound cannot be viewed as being 
absolute. Nevertheless, we believe it can provide a use¬ 
ful range of scales for experimentalists which makes it 
worthwhile for further investigation. 


VI. MINIMAL LEPTOBARYONS AT THE LHC 

To definitively test this model, one would ideally like 
to produce the new leptobaryons 0, E, y at a collider and 
measure their electroweak quantum numbers. However, 
the requirement imposed by dark matter relic abundance 
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FIG. 7. Upper limits on an as a function of Zb mass drawn 
from the absence of a dijet resonance signal. Note that ob is 
on a logarithmic scale. 


discussed in the previous section suggests that the light¬ 
est leptobaryon y is more likely to have a mass in the TeV 
range. On account of the accidental global Zi symme¬ 
try, leptobaryons must be produced in pairs. This leads 
to phase space suppression in the production rates, and 
makes their direct observation more difficult. 

However, it is possible to reveal the underlying bary- 
onic symmetry breaking mechanism by producing and 
studying the gauge boson Zb and the Higgs boson 5'b, 
which will be discussed in this section. Although these 
new bosons are not specihc to this model, for very small 
values of the Higgs mixing angle |0| < 10“^, we find that 
it may be possible to indirectly observe the charged lepto¬ 
baryons based on their imprints left on the loop-mediated 
decay modes of S'b. 


A. Production of Zb and limits on ob 


quark pairs as the dominant mode^ The gauge boson 
Zb universally couples to all flavors of quarks with equal 
strength, implying equality among the branching frac¬ 
tions to quark pairs. Therefore, the appropriate search 
strategy at the LHC is to look for a dijet formed by the 
lighter quarks (u, d, c, s, b) in the final state, along with 
a corresponding absence of a dilepton resonance to verify 
the leptophobic nature of Zb. 

Assuming massless quarks, the leading order parton- 
level cross section for a given flavor of initial and final 
state quarks is 


dTrc^B s 

243 {s - + ml^Tzs 


mz^ , 
243 Tz^ ^ 




(nwa) 


(42) 


where is the partonic center-of-mass energy and 
= 2aB'mzjil^ is the Zb total width. In the sec¬ 
ond line, the narrow width approximation is made and 
is justified because <C mzg |23]. The partonic 

cross section is convoluted with the leading order MSTW 
2008 [30] parton distribution functions for -y/s = 8 TeV, 
and summed over the five light quark flavors yielding an 
estimate of the pp jj cross section. 

The most stringent bound on the acceptance times 
cross section in the range 1 TeV < m^g < 3 TeV comes 
from the recent CMS search for a dijet resonance |3T] with 
19.7 fb“^ of data. Using their suggested acceptance fac¬ 
tor of 0.6, we convert their limits into upper limits on ctb 
as a function of Zb mass, shown in Fig. Also shown in 
dashed lines are the more competitive limits in the lowest 
mass region [32] derived from the CDF (0.11 fb~^) [33] 
data set and an older CMS analysis (4 fb“^) |5i] . 

We find that constraints on the coupling constant de¬ 
pend strongly on the Zb mass. Very roughly, the cou¬ 
pling is constrained to lie below ub ^ 0.02, except for 
the narrow mass range 1.2 TeV < m^g ^1.7 TeV where 
the upper bound is at an ~ 0.008. The dark matter pa¬ 
rameter scan of the previous section is indicated by the 
dashed green rectangle, and lies just below the ub upper 
limit set by the recent CMS analysis. 


We identify two dominant production modes of Zb at 
the LHC: the s-channel production pp jj leading 
to a dijet resonance, and the production in association 
with the baryonic Higgs pp Zb 5'b. Here we discuss 
the resonant Zb production with which we place a limit 
on the mass m^g and coupling q;b based on the absence 
of a Zb signal at the LHC, and leave the discussion of 
associated production to the next subsection. 

Because dark matter relic abundance considerations 
suggest > mz^l^ (see Fig. [^, the decay channels 
to leptobaryon pairs 

Zb ^'tpilj,T,T,,xx (41) 

are closed. This leaves the decay into standard model 


B. Decay and production of baryonic Higgs Sb 

As long as they are open, the largest decay modes of 
the baryonic Higgs 5 b for small 9 are to a pair of lepto¬ 
baryons or to a Zb pair 

5 b ijij, SS, XX, Zb Zb . 


^ Subdominant modes induced by fermion loops to final state Higgs 
bosons H and 5 b, and to electroweak gauge bosons W, Z, and 7 
are possible. A subset of these are studied in m in the context 
of indirect detection in model VA. 
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FIG. 8 . Representative Feynman graphs contributing to 
baryonic Higgs decay to electroweak gauge bosons V, V' = 
{W^, Z, 7 }. The tree-level graph on the left is mixing angle 
suppressed compared to the loop diagram on the right. 


Since the requirement of dark matter relic abundance 
saturation favors heavy dark matter in the non-resonant 
regime (m^ > tosb/ 2 ), we expect the leptobaryon decay 
modes to be closed. This leaves the Zb pair as the single 
dominant decay channel for small mixing angles 6 < 0 . 1 . 
For larger mixing angles, decay modes to the ZZ, WW 
and HH bosons are dominant. 

Below the ZbZb threshold, the baryonic Higgs decay 
profile exhibits a striking interplay between the mixing 
angle suppressed and the loop suppressed decay modes. 
Particularly interesting are decays to electroweak gauge 
bosons shown in Fig. The dramatic behavior of the 
branching fractions of 5 'b as a function of the mixing an¬ 
gle 9 is illustrated in Fig. [^and is fairly independent of 
the particular values taken for the masses ms, mzj^, m^, 
coupling Ob and sign of the mixing angle. For simplicity, 
masses of the charged leptobaryons have been taken to 
infinity (soft Sb limit); we discuss the impact of finite lep¬ 
tobaryon masses consistent with vacuum stability shortly 
below. 

For mixing angles moderate in magnitude |0| > 10“^, 
in addition to the HH mode, the baryonic Higgs in¬ 
herits the largest decays of the standard model Higgs, 
which are the tree-level modes WW, ZZ and tt. As the 
mixing angle is parametrically decreased, the standard- 
model-like decay modes give way to the leptobaryon loop- 
mediated decay modes. In the small mixing angle regime 
\6\ < 10 “^, the modes to electroweak gauge bosons dom¬ 
inate the total width. The loop-mediated decay to a top 
quark pair is suppressed by due to chiral sym¬ 

metry. 

The relative strengths of the decays to vector bosons in 
the limit 6 * —>■ 0 are easily understood by integrating out 
the charged leptobaryonic multiplets dt and S. The effec¬ 
tive Lagrangian describing the coupling of the baryonic 
Higgs field s to the electroweak gauge bosons is 


^eff = 


gB 


5] [e^Q) 

f 


+ ^{Tj-QfslfsZ^,,Z^'' 

Cw 

+ ^Qf{Tf-Qfsl)sF^,Zf^'' 

Cw 


9b 






(43) 


where the sum in the first three terms is over the isospin 


ag=0.005, m^ = l TeV, 

=1.2 TeV, = 1.2 TeV, m-^, m,p —> 00 



FIG. 9. Branching fraction of the baryonic Higgs boson Sb 
as a function of the mixing angle 6 \6 > Q solid, 9 < 0 dashed] 
below the 5 'b —>■ Zb Zb threshold. For mixing angles moderate 
in magnitude \6\ > 10“® standard model Higgs-like decay 
modes dominate, and for smaller mixing angles \6\ < 10 “^ 
loop-mediated modes dominate (see Fig. |^. The vertical line 
at \6\ = 0.35 is the upper limit determined from Higgs signal 
strength measurements. 


components / = "0°, S’*", S°}, and the sum in the 

final term runs over entire multiplets F = {'F, S} where 




=1 and = 2. The relative partial widths are 
obtained by squaring the couplings, and dividing the ZZ 
and 77 modes by 2 to account for indistinguishability of 
final state particles. This yields the relative branching 
ratios of the baryonic Higgs [see (C2)] 


T 


WW 


'Wzz ■ F^^ : r.vT' = 20 : 7 : 3 : 1. 


(44) 


It is particularly interesting to compare this prediction 
with that of leptobaryon model VA. The predicted ratios 
derived from the electroweak currents in Eq. (Al) are 


F 


WW 


■rzz-r 


Zj 


: = 2 : 1 : 10 


-3 


: 1. (model VA) 


The qualitative profile is substantially different from 
that of model A, with the Zj mode being far below the 
77 mode in model VA. Note that the precise numerical 
ratios are subject to mass corrections. The largest 
corrections come from taking finite values of the charged 
leptobaryon masses to be consistent with the vacuum 
stability bound in Section |HIB| However, the gross 
hierarchy of decay modes remains largely unchanged. 
Therefore, for small mixing angles, we are optimistic 
about the prospects of distinguishing the leptobaryon 
models, even if the fermions cannot be directly produced 
at the LHC. 


We turn to the production of the baryonic Higgs at 
the LHC. Similar to its decay, the production cross sec¬ 
tions depend on the mixing angle 9 and also exhibit an 
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a'g=0.005, m^ = l TeV, M^^=1.2TeV 



ff5=0.005, 6»=±0.05, Mz^ =1.2 TeV 



mg [GeV] 


FIG. 10. Baryonic Higgs production cross sections at a gfs = 14 TeV pp collider as a function of mixing angle (left) and Sb 
mass (right) [6 > 0 solid, 6 < 0 dashed]. The curve labeled “VBF” is the combined contributions of the mass suppressed ZbZb, 
the mixing angle suppressed WW/ZZ, and the loop suppressed WW, ZZ, Zy, and 77 f usio n (see Table III I. The effect of 
charged leptobaryon loops is included by inserting the appropriate effective operator in Eq. (431. The vertical line at | 6 '| = 0.35 
is the upper limit determined from Higgs signal strength measurements. 


interplay between mixing angle and loop suppressed con¬ 
tributions. The production channels and representative 
parton-level Feynman graphs are listed in Table Ell 


Using the effective Lagrangian in Eq. (43) in the limit 
—>■ 00 , we included the loop-mediated leptobaryon 
contributions to the vector boson fusion (VBF) and as¬ 
sociated production channels. In order to correctly han¬ 
dle interference and to facilitate the phase space integra¬ 
tion, we obtained all production cross sections, except for 
gluon fusion, with CalcHEP 3.4 [3S], using the leading 
order MSTW 2002 parton distributions. To keep finite 
top quark mass dependence, the one-loop gluon fusion 
cross section was computed separately with Package-X 
[36] . and the results were convoluted with the MSTW 
2008 PDF. We checked our CalcHEP implementation 
of leptobaryon model A by verifying that the tree-level 
predictions for and ZSb production cross sections 

computed by hand and convoluted with the MSTW 2008 
PDFs matched the CalcHEP results. This also con¬ 
firmed that the results are fairly independent of the PDF 
data sets used. 


The behavior of the cross sections is shown in Fig. 
for fixed ae = 0.005 and ttizb = 1-2 TeV. The panel on 
the left illustrates the mixing angle dependence of the 
various cross sections at fixed mg = 1 TeV. At small 
mixing angles, the tree-level ZbSb cross section is the 
largest at 0.05 fb. But, for larger mixing angles, the 
VBF and gluon fusion modes become larger. The panel 
on the right shows the mg dependence at the mixing 
angle |0| = 0.05. For light baryonic Higgs, the dominant 
mode of production is gluon fusion with a cross section 
reaching 20 fb. 

The LHC search strategy for the baryonic Higgs is in¬ 
fluenced by the production cross section and decay pro¬ 
file. For mixing angles below 10“^, the produc- 


Associated production 



Vector boson fusion (VBF) 

qq —> qqSB, qq —> qqSB, qq —qqSB 



gg —Sb 



TABLE III. Production channels of the baryonic Higgs Sb at 
a hadron collider. Vertices marked with a ‘0’ are suppressed 
by the mixing angle, and square vertices are derived from the 
effective Lagrangian in Eq. (43l. 


tion cross sections are quite small (0.05 fb), but the rel¬ 
atively clean 77 , Zy and ZZ decay modes are large and 
may make it possible to identify a signal at the LHC by 
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searching for a resonance in the appropriate channel. The 
current ATLAS limit on crfld(*S'B) • BR(S'b —t 77 ) [37] at 
•y/s = 8 TeV for the mass range 100-600 GeV is around 
5 fb (95% C.L.), and is currently too large to place mean¬ 
ingful bounds. For larger mixing angles |0| > 0.02 the 
gluon fusion mode is very large, although the only rela¬ 
tively clean decay mode is the ZZ. 

However, by LHC run 3, the CMS and ATLAS exper¬ 
iments are expected to collect 1000-3000 fb“^ of data at 
y/s = 14 TeV. This will enable the experiments to probe 
the small mixing angle scenario for masses of Sb and Zb 
reaching around 1 TeV (Fig. 10 left panel). Therefore, 
we expect that the experiments at the LHC will be able 
to probe this model in the near future. 


VII. CONCLUSIONS 

In this article, we initiated an investigation of the dark 
matter and LHC phenomenology of low scale U(l)g sym¬ 
metry breaking of leptobaryon model A. Anomaly can¬ 
cellation demands the existence of new colorless fermions 
■ 0 , S and X carrying baryon and lepton numbers, called 
leptobaryons. The lightest leptobaryon y is a Majorana 
fermion, and is automatically stabilized by an acciden¬ 
tal Z 2 symmetry making it a dark matter candidate. In 
this model, saturation of the observed dark matter relic 
abundance favors heavy dark matter (m^, > 1 TeV), and 
null results by direct detection experiments lead to an 
upper limit on the Higgs mixing angle (| 6 *| < 0.22) that 
is competitive with current limits from LHC experiments 
(1^1 ^ 0.35). Future experiments will tighten these limits. 
One-loop vacuum metastability places an upper limit on 
the leptobaryon masses not far above the baryonic gauge 
and Higgs boson masses. Although it may be difficult 
to produce and study the leptobaryons at the LHC, it is 
possible to discover the U(l )3 symmetry breaking mecha¬ 
nism by producing the baryonic gauge and Higgs bosons. 
For small mixing angles, the branching fractions of the 
leptobaryon loop-mediated ^b —t 77 , ^ 7 , ZZ modes can 
reveal the underlying fermion content at the LHC. 

Our initial study opens up new avenues for future in¬ 
vestigations. It may be interesting to speculate whether 
the other neutral leptobaryons 0 ° or can populate 
the dark matter abundance. Furthermore, in the limit 
of vanishing Higgs Yukawa couplings Ai... A 4 in Eq. ^ 
all neutral leptobaryons are stable, and give rise to mul¬ 
ticomponent dark matter. Within the context of elec- 
troweak baryogenesis, the addition of the baryonic Higgs 
leads to another scalar field participating in the elec- 
troweak phase transition. Furthermore, the standard 
model Higgs coupling to the leptobaryons can accommo¬ 
date new sources oi CP violation. It may be worthwhile 
to investigate whether this model can account for the 
observed baryon asymmetry while evading constraints 
from low energy electric dipole moment measurements. 
Finally, spontaneous symmetry breaking of U(l) gauge 
groups leads to the generation of Nielsen-Olson strings 


Gauge fields 

SU(2)^ 

we 

U(l)y 

Re 

U(1)b 

U(1)l 

Fermions 


1 

0 

0 

-1 

T 

2 

-1/2 

Hi 

Li 

T 

2 

1/2 

-B2 

-L2 

V 

1 

-1 

B2 

L2 

V 

1 

-1 

-Bi 

-Li 

X 

1 

0 

B2 

L2 

X 

1 

0 

-Hi 

-Li 

Scalar fields 

H 

2 

1/2 

0 

0 

5'b 

1 

0 

-3 

-3 

5l 

1 

0 

0 

2 


TABLE IV. Field content of leptobaryon model VA, with the 
fermions given in two-component left-handed notation. All 
fields are color singlets, and the hypercharge is normalized 
according to Q — + Y. 


in the early Universe. In light of the recent study by 
[551155] , it may also be interesting to investigate whether 
high scale U( 1 )l breaking can lead to observable signa¬ 
tures. 
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Appendix A: Leptobaryon model VA 


For reference, the formulation of leptobaryon model 
VA is briefly summarized and some salient phenomeno¬ 
logical features of low scale U( 1)0 breaking are high¬ 
lighted in this appendix. 

The gauge-Higgs sector is identical to that of model A. 
But, in model VA six fermionic leptobaryon multiplets 
are added (see Table [rV]) . Upon spontaneous symmetry 
breaking of U( 1 ) 0 , the leptobaryons acquire vectorlike 
masses leading to three Dirac fermions: (0'^, 0~), r]~ and 
X- They couple to the standard model through the weak 
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and electromagnetic vector currents given by 

•^NC = jV'W'” + (-5 + + slr]-J^^r]- 

Jcc = . (Ai) 

and through the extended baryonic and leptonic currents, 
which have both vector and axial-vector parts 


\ 7M 

'^B,ext — V “t" '^R 


B,V 


^B,A ’ 


tM _ tM 

^L,ext ~ ^L,V 


■ 


(A2) 


The axial-vector parts are fixed by the anomaly cancel¬ 
lation condition {Bi — B 2 = Li — L 2 = —3): 


•^B.A = Jl,A = I (V'°7^75V'° + ) 

- §7^7^757” - ix7''75X • (A3) 

But the overall normalization of the vector parts remains 
unconstrained 

■^B,v = + xYx) 

Jl V = + xi^x) , 

where B = Bi + B 2 and L = Li + L 2 . 

The fact that the leptobaryons have a vector part to 
their interactions with Zb qualitatively changes the dark 
matter phenomenology as compared to that of model A. 
The direct detection cross section via t-channel Zb ex¬ 
change is not velocity suppressed. Also, the resonant 
annihilation channel XX QQ is not velocity sup¬ 

pressed. This makes it possible for lighter dark matter, 
and opens the possibility for signals in indirect detection 
experiments [28) . Furthermore, the vector currents are 
associated with an accidental global “77 symmetry” that 
is anomaly free [43], and leads to an interesting connec¬ 
tion between the observed baryon asymmetry and the 
dark matter asymmetry [7|. 


Appendix B: Scalar potential parameter relations 

For reference, we provide the tree-level relations be¬ 
tween the scalar potential parameters in 0 and the 
masses, VEVs and mixing angle 

cos^ 9 + to| sin^ 6) 

- -ml)s\TL{29) 

sin^ 0 + rrig cos^ 0) 

A = cos^ 0 + m| sin^ 0) 

b = sin^ 0 + rrig cos^ 0) 

« = ) sin(20). (BI) 


Appendix C: Baryonic gauge and Higgs boson 
partial widths 


The partial widths of the baryonic gauge boson are 
given by 


r(ZB 


aB 

m) = 

3aB 



1-4- 


m 


Zb 




3/2 


(Cl) 


To avoid displaying the complicated interference 
terms in the 5 'b partial widths, we give the formulas in 
two distinct limits: 0 —)■ 7r/2 and d —)■ 0. In the limit 
0 —7r/2, the tree-level standard model Higgs-like modes 
are dominating 

r(S'B ^ WW) = 

Sv^TT 
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In the limit 9 —r 0 the only tree-level mode is ^b —>■ 
Zb Zb. The leptobaryon loop-mediated modes to elec- 
troweak gauge bosons are given in the limit rnw,z,'y ‘C 
ms 'C rn^ Y., and can be derived from Eq. (43) 


r(5'B —>■ ZbZb) = 
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The decay mode to HH vanishes in both limits. Here 
we give the tree-level partial width that is valid for all 
values of 0 


r(^B ^ HU) = 


where 


CSHH = 


[2m\ + ml) 



(C3) 

oTTUis y Trig 

■ /„„x/cos 0 sin 0 \ 
sin 26» --h- . 

\ V Vb / 

(C4) 


These formulas are valid for \9\ < ir/i. 
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